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AQVAL/l  (AQ7) 
User's  Guide  and  Program  Description 


ABSTRACT 


AQVAL/l  (AQ7)  is  a  PL/l  program  which  synthesizes  quasi- 
optimal  formulas  of  the  variable-valued  logic  system  VL-.  .   By  'quasi- 
optimal  formulas'  we  mean  here  disjunctive  simple  VL]_  formulas,!  which 
are  optimal  or  sub-optimal  with  regard  to  a  user- specified  optimality 
functional. 

The  basic  application  of  the  program  is  in  the  area  of  machine 
learning  and  inductive  inference  ('inductive  learning'):   from 
descriptions  of  objects  with  known  class  membership,  the  program  infers 
optimal  or  sub-optimal  descriptions  of  object  classes.   These  descriptions 
are  expressed  as  VL]_  formulas  and  represent  certain  generalizations  of 
inputted  information.   The  program  can  also  be  used  (at  the  appropriate 
setting  of  its  parameters)  for  an  efficient  minimization  of  binary-  or 
multi -valued  switching  functions  with  a  large  number  of  variables  (e.g., 
50-100  variables) . 
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AQVAL/l:   Part  I 
1.   INTRODUCTION 

This  paper  provides  a  user's  guide  and  brief  description 
of  the  program  AQVAL/l  version  AQJ.   The  program  infers  an  optimal  or 
sub-optimal  (with  regard  to  a  user- specified  optimality  functional) 
disjunctive  simple  formula (s)  of  the  variable -valued  logic  system  VL  . 

y  y 

In  this  version,    the  formulas  are  inferred  from  event  sets  (i.e., 
sets  of  sequences  of  values  of  input  variables)  where  each  set  is 
associated  with  a  certain  value  of  the  output  variable(s).   The  way 
in  which  sets  are  associated  with  values  of  output  variable (s)  depends 
on  the  mode  of  program  operation,  as  explained  in  section  2.2.  in  the 
description  of  parameter  MODE. 

It  is  assumed  here  that  the  reader  is  acquainted  with  the 
definition  of  the  variable-valued  logic  system  VL  given  in  paper1  and 
therefore  very  limited  explanation  of  the  concepts  of  VL  is  included. 

Chapters  2  and  3  describe  in  detail  the  input  and  output 
parameters  of  the  program.   Chapter  k   gives  two  complete  examples  of 
actual  specification  of  these  parameters.   Chapter  5  describes  briefly 
the  implementation  of  the  synthesis  algorithms  (described  in  paper2) .  in 

order  to  understand  the  last  chapter,  the  reader  should  be  acquainted 

o 

with  paper  .   Bibliography,  in  addition  to  papers  containing  the  background 

information  for  this  paper  (i.e.,  the  definition  of  the  VL  and  a 
description  of  synthesis  algorithms2)  includes  also  the  complete  list  of 
papers  relevant  to  the  current  status  of  variable-valued  logic  and  its 


*-* 


The  name  AQVAL/l  was  derived  from  'Algorithm  A-  applied  for  the  synthesis 
of  Variable-Valued  Logic  formulas. 

A  synthesis  of  quasi-optimal  formulas  from  other  (non-optimal)  formulas 
(rather  than  from  event  sets)  is  presently  implemented  only  in  version 
AQ9  of  AQVAL/l  described  in  paper20. 


applications .   In  Appendix  A  a  short  form  describing  steps  in  input 
preparation  is  given  together  with  a  list  of  optimality  criteria. 

The  described  here  version  AQ7  is  written  in  PL/l  language* 
and  consists  of  approximately  860  statements.  The  object  module 
requires  about  60k  bytes  (together  with  the  PL/l  library  routines 
about  90k)  plus  storage  for  the  data,  which  may  take  somewhere  between 
10-150k  bytes,  depending  on  the  total  number  of  events,  number  of 
variables,  sizes  of  the  domains  of  variables,  and  some  other  parameters 
(minimum  ~112k  bytes) .  The  complete  listing  of  the  program  is  given  in 
Appendix  B. 

2.   INPUT  PARAMETERS  AND  FORMAT 
2.1  Description  of  input  parameters 

Following  is  a  description  of  all  input  parameters  in  the 
order  and  format  in  which  they  should  be  specified  in  the  input  stream 
to  the  program.   The  description  of  each  parameter  consists  of: 
i.   the  name  of  the  parameter, 
ii.   an  example  of  its  specification 
iii .   a  list  of  possible  values  of  the  parameter 
iv.   a  description  of  its  meaning. 
Input  parameters  are  divided  into  2  groups : 
A.   Control  parameters  (specified  in  the  PL/l  input  data  format)  which 
supply  an  information  about  the  way  the  program  should  be  executed. 

Specification  of  these  parameters  ends  with  a  semi-colon  (;). 
If  any  control  parameter  is  not  specified  by  the  user  in  his  input 
data  then  it  takes  a  default  value.   The  default  value  is  indicated 
as  the  value  used  in  the  example  of  the  parameter  specification. 


It  is  adapted  for  optimizing  compiler  and  runs  under  the  HASP 
operating  system. 


B.   Data  parameters  which  describe  the  VL  function(s)  to  be  optimized 
and  optimization  functional. 

2.2  Steps  in  the  input  specification: 

I.   Specification  of  control  parameters: 
•  MODE 
Example : 


MODE  =  'IC1 


Possible  values:   'IC,  'DC,  'VL'. 

This  parameter  defines  the  mode  of  program  operation 
that  is,  it  defines  how  the  program  should  interpret  the 
data  parameters  (specifically,  event  sets).   Consequently, 
the  interpretation  of  the  output  from  the  program  will 
also  be  dependent  on  this  parameter. 

Here  we  will  present  a  brief  description  of  each  mode 
(IC,  DC,  VL),  and  in  Chapter  k-   we  will  illustrate  each 
mode  by  an  example.   (For  a  more  detailed  description  of 
these  concepts,  consult  paper  . ) 

For  the  convenience  of  the  description  of  MODE  parameter, 
let  us  assume  that  in  the  input  stream  there  are  m  (disjoint) 


event  sets  (step  TV,  p.  Ik,    parameter  F)  which  we  will  denote 
here  as  fO,  Fl,  .  . . ,  Fm" -1- .   (The  assignment  of  symbols 
F°,  F1,  . .  .,Fm~l  to  the  event  sets  specified  in  the  input 
stream  depends  on  the  OLIST  parameter  (step  TV,  p.  ll)). 


--  If  MODE  =  'IC  ('Intersecting  Covers'),  then  the 
event  sets  are  treated  as  defining  a  set  of  m 
binary-valued  VL  functions  fj_,  i=0, 1,  .  .  .,m-l. 
A  function  fj_  takes  value  1  for  events  from  F^, 
and  value  0  for  the  remaining  events,  specified 
in  FJ,  jVi,  j=0,l,  ...,m-l: 


where  F  Is  the  union  of  all  events  specified 
in  the  program,  i.e. 

F  =  F°  U  F1  U  . . .  U  F*-1 

The  output  from  the  program  consists  of  m  DVL^ 
formulas  which  are  optimized  expressions  V(fi) 
(with  regard  to  the  optimality  functional  A, 
step  IV,  p.  11,  parameter  NCRIT  CLIST(NCRIT) 
TLIST(NCRIT))  of  functions  f . . 


These  expressions  generalize  the  input 
information  about  functions  in  the  sense  that 
they  may  now  assign  a  specific  value  (l  or  0) 
to  events  which  were  not  listed  in  the  input 
stream  (i.e.,  events  from  e\f,  called  *-events, 
where  E  is  the  universe  of  events) .   Formally, 
a  V(fj_)  expression  is  defined  as  a  DVL]_  expression 
corresponding  to  an  optimized  (with  regard  to 
functional  A)  cover  of  set  F*  against  the  set 
F\Fi: 

V(fi)  *-C±  =   C(F1|F\Fi)   1=0,1,  ...,m-l 

In  this  mode  some  ^-events  can  satisfy  more 
than  one  formula  V(fj_),  i.e.,  the  intersection  of 
some  formulas  V(fj_)  may  not  be  empty  (hence  the 
mode  is  called  'Intersecting  Covers').   Note, 
however,  that  each  event  from  F(i.e.,  event  from 
the  input  list)  satisfies  exactly  one  formula 
(it  follows  from  the  assumption  that  sets  F3-  are 
disjoint) . 

If  MODE  =  'DC  ('Disjoint  Covers'),  then  the  input 
event  sets  are  treated  similarly  as  in  the  case 
when  MODE  =  'IC.   The  difference  is,  however,  that 
the  program  will  produce  formulas  V(fj_)  which  are 
now  all  disjoint  (a  given  *-event  will  satisfy  at 
most  one  output  formula  V(fi)).   Formulas  V(fj_), 
i=0, 1,  .  .  .,m-l,  obtained  in  this  mode,  correspond 
to  consecutively  synthesized  covers  C^ : 

m-1   . 
V(fQ)   —   C0   =   C(F°  |  (J  FJ) 


-,        m-1 
V(f1)    —   C1   =   C(FL  |  CUQ   U  (J  FJ) 

3=2 


i-1      m-1 
V(f.)    4-   C    =   C(FX  |  (J  CV  U   (J   FJ) 

j=0  J   j=i+l 


..   m-2 

j=0  J 


(By  C.  is  meant  the  union  of  complexes  in  cover 
Cj.)  J 

Each  formula  V(fj_)  is  a  union  of  terms,  each 
term  corresponding  to  a  complex  in  the  cover  Cj_ . 

As  we  see,  the  formulas  depend  on  the  assumed 
order  of  sets  Fi  (see  parameter  OLIST  in  step  IV 
p.  11). 

If  MODE  =  'VL',  then  sets  F1,  i=0,l,2,  . .  .,m-l  are 
treated  as  sets  defining  one  m-valued  VL  function: 

f :  E  ->  {0,1,  .  ..,m-l} 

where 

E  --  the  event  space  defined  as  D^  XD2  x  .  ..  xD  , 
(Dj_  -  domains  of  individual  variables) 

Sets  F1  are  related  to  the  function  as  follows: 

F1  =  {e  eE|f(e)  =  i)   1=0,1,  . .  .,m-l 

The  output  from  the  program  is  an  optimized  DVL-, 
formula  V(f)    expressing   function  f.      Formula  V(f) 
is  defined  as: 

v(f)  =  MK^v  i^v  ...)  v 

(m-2)(T^"2  V  T^'2   V  .-.)  V 


1  (t\   V  T^  V  ... 


where 
1-] 


T7  ,  i=l, 2, ...  --  terms  corresponding  to  complexes 

in  the  quasi-optimal  (q-o)  cover 
m-2 

c_tl  =  cCf"1"1!  u  f±) 


m-  _ 

i=0 


T.   ,  1=1,2,...  --  terms  corresponding  to  complexes  in 

the  q-o  cover 

C    =   C(Fm"^|  U  F.) 

m"2        1=0  x 


T.,  i=l,2, ...   --  terms  corresponding  to  complexes  in 


1 


the  q-o  cover 

c1  =  c(f1|f°) 


Note  that  the  total  number  of  covers  C.  is  smaller  by 
1  than  the  number  of  sets  F1 • 


•   INFORM 


Example 


INFORM  =  'VECTOR' 


Possible  values:   'VECTOR',  'GAMMA',  'VL'* 

If  INFORM  equals  'VECTOR',  then  events  are  read  in  as 
sequences  of  variable  values.   If  INFORM  equals  'GAMMA.',  then 
events  are  read  in  the  form  of  event  numbers  7(e) .   If 
INFORM  =  'VL',  then  events  are  expressed  as  a  VLn  formula. 


•  TITLE 
Example : 


TITLE  =  3 


Possible  values:   any  non-negative  integer 

This  parameter  specifies  the  number  of  lines  reserved  for 
the  problem  title  (which  is  typed  just  after  the  semi-colon 
which  ends  the  specification  of  control  parameters). 


MAXSTAR,  CUTS  TAR 


Example : 


MAXSTAR  =  150 
CUTSTAR  =  50 


Possible  values 


MAXSTAR  I 
CUTSTARJ 


any  positive  integer  smaller 

than  NGE  value  (see  next  parameter) 


These  parameters  are  used  to  control  the  speed  of  the 
program.   Small  values  of  these  parameters  (e.g.,  MAXSTAR  =  5, 
CUTSTAR  =  3)  will  cause  the  program  to  run  more  quickly  (but 
may  also  cause  the  program  to  generate  less  optimal  formula (s)) 

The  work  of  these  parameters  can  briefly  be  described  as 
follows:  A  star  generation  is  a  multi-step  process  which  at 
each  step  produces  a  set  of  complexes  ('intermediate  star'). 
If  at  any  step  an  intermediate  star  contains  more  than  MAXSTAR 
complexes,  then  it  is  reduced  to  CUTSTAR  complexes,  selected  as 
the  'best'  according  to  the  assumed  optimality  functional  (see 
CLIST  and  TLIST  parameters) . 


■  op) 

Implemented  presently  only  in  version  AQ9  described  in  paper^  . 


NGE 

Example:    NGE  =  200 

Possible  values:   any  positive  integer  (smaller  than  2  ^) 

This  parameter  denotes  the  maximal  number  of  complexes 

which  can  be  stored  during  a  star  generation.   NGE  defines* 

the  maximum  size  of  a  list  G  each  element  of  which  consists 

of  k   fields: 

n 

1.  storing  a  complex  (which  is  a  bit  string  of  length  Z  d. ) 

i=l  x 

2.  storing  a  value  (real  number)  of  a  criterion 

under  consideration  for  the  complex 

3 •  and  k .   pointer  to  the  next  and  to  the  previous  element 
of  the  list. 


•  NMQ 


Example : 


NMQ  =  25 


Possible  values:   any  positive  integers 

This  parameter  defines  the  length  of  a  list  reserved 
for  storing  complexes  L^  ( 'quasi-extremals '  or  'best' 
complexes)  constituting  the  final  formula(s). 


•  LQST 


Example : 


LQST  =  '1 


Possible  values:   '0*B,  '1'B 

If  LQST  =  ' 1'B  each  complex  L^  selected  from  a  star 
will  be  'reduced'  before  it  is  put  on  the  output  list. 
The  reduction  means  that  the  references  (ranges)  of 
variables  in  each  complex  will  be  reduced  as  much  as 
possible,  providing  that 

1.  the  complex  will  still  cover  the  same  number  of 
events  in  the  original  set  F1 

2.  the  reduced  complex  will  have  the  same  number  of 
literals  (i.e.,  selectors  in  the  corresponding 
term  is  the  formula) . 


See  page '8  of  the  program  listing  in  Appendix  (separate  report) 


SAVE 

Example:     SAVE  =  '0*B 

Possible  values:   'O'B,  'l'B 

If  SAVE  =  'O'B,  the  output  formula  is  not  stored. 
If  SAVE  =  'l'B,  then  the  formula  is  stored  in  the  disc  file 
called  'COVER'  (the  disc  file  is  specified  by  JCL  cards 
listed  with  example  1  at  the  end  of  this  paper) . 

Individual  complexes  are  stored  as  bit  strings 
('binary-positional'  notation)  in  PL/l  list  format  (see  Chapter  5), 
similarly  as  they  are  stored  during  the  program  execution. 
This  way  of  storing  facilitates  the  evaluation  of  the  formula. 


•  QLQT 


Example : 


QLQT  =  'O'B 


Possible  values:   'O'B,  'l'B 

This  parameter  controls  the  printing  of  summary 
information  about  star  generation.   If  QLQT  equals 
'O'B,  no  information  will  be  printed.   If  QLQT  equals 
'l'B,  then  the  program  will  print  out  a  summary  about  each 
star,  the  qua si -extremal,  and  the  number  of  new  events  covered 
by  it  in  a  given  F1 . 


LQTEACE 


Example 


LQTRACE  =  'O'B 


Possible  values:   'O'B,  'l'B 

This  parameter  controls  whether  detailed  information 
about  each  quasi -extremal  should  be  printed.   If 
LQTRACE  =  'O'B,  then  no  information  is  printed.   If 
LQTRACE  =  'l'B,  then  in  addition  to  the  information 
printed  by  QLQT,  the  actual  cost  of  the  quasi-extremal 
and  the  array  denoting  the  elements  not  covered  by  the 
current  and  previous  stars  is  printed. 


•  STRACE 


Example 


STRACE  =  'O'B 


Possible  values:   'O'B,  'l'B 

This  parameter  controls  whether  detailed  information 
about  star  generation  should  be  printed.   If  STRACE  =  'O'B, 
then  detailed  information  is  not  printed.   If  STRACE  =  'l'B 
each  step  of  star  generation  is  printed. 


•  QST 

Example : 


QST  =  'O'B 


Possible  values:   'O'B,  'l'B 

This  parameter  controls  whether  a  summary  of  the  steps 

in  star  generation  should  be  printed.   If  QST  =  'O'B,  then 

the  summary  is  not  printed.  If  QST  =  'l'B,  then  the  summary 
is  printed. 


II.   Separation  between  control  and  data  parameters: 

Type  semi-colon  (;)  after  the  last  control  parameter.   If  no 
control  parameters  are  specified  (i.e.,  their  default  values  are 
accepted),  the  semi-colon  still  should  be  typed. 


III.   Specification  of  the  problem  title: 

Type  the  title  of  the  program  using  as  many  lines  as  were 
specified  by  TITLE  control  parameter.   The  title  information  will 
be  reproduced  verbatim  at  the  top  of  the  output  page. 


TV.   Specification  of  data  parameters 

All  of  the  following  data  parameters  are  in  PL/1  LIST 
format;  that  is,  a  value  must  be  present  for  each  of  these 
parameters  in  the  order  specified.   There  are  no  default  values 
for  these  parameters . 


NSPEC  TYPE  TYPELIST ( NSPEC' 


Example 


2  'FACTOR'  1  k 


Possible  values 


NSPEC  —  a  positive  integer  g  6k    (this 

restriction  can  be  easily  removed) 

TYPE  --  'FACTOR'  or  'INTERVAL' 

TYPELIST  --  a  list  of  NSPEC  positive  integers  g  6k 

These  parameters  specify  the  type  of  each  variable  which 
is  used  in  the  event  description.   NSPEC  specifies  the  number 
of  variables  which  will  be  of  the  type  given  by  the  parameter 
TYPE.   The  parameter  TYPE  specifies  the  type  of  the  variables 
whose  indices  are  specified  by  TYPELIST.   (All  variables  not 
in  TYPELIST  are  assumed  to  have  the  'opposite'  type.) 
There  should  be  NSPEC  indices  in  TYPELIST.   In  the  example 
above,  there  are  NSPEC  =  2  variables  which  are  given  the  type 
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'FACTOR'   (i.e.,  xl  and  xK) .   The  rest  are  given  the  'opposite' 


type  namely  'INTERVAL*   (i.e.,  x2,  x3,  x5, 
variable  indices  begin  with  1. 


. ) .   Note  that 


The  following  are  some  guidelines  for  selecting  the  type 
of  each  variable. 

If  a  variable  takes  values  which  have  natural  linear 
order  (e.g.,  represent  temperature,  height,  length,  grey- 
level  of  a  picture  element,  etc.)  then  it  is  advisable  to 
specify  it  as  an  interval  variable  ('INTERVALS'  type).   If  this 
is  not  the  case  (e.i 


■  •  > 


the  variable  represents  a  set  of 


independent  objects,  relations,  properties,  etc.,  i.e.  is 
measured  on  a  nominal  scale),  then  the  variable  should  be 
specified  as  a  factor  variable  ('FACTORS'  type). 

It  should  be  noted  that  interval  variables  are  simpler 
for  evaluation  than  factor  variables.  Also,  the  synthesis 
of  optimal  VL]_  formulas  usually  will  take  considerably  less 
computational  time  and  memory  if  variables  are  specified  as 
interval  rather  than  as  factor  variables .   On  the  other 
hand,  VLi  formulas  with  interval  variables  will  usually  have 
considerably  more  terms  and  selectors  than  formulas  with 
factor  variables  (when  no  restrictions  on  the  variables  are 
assumed) . 


NV  NLEV  NV. 


Example : 


^  5,  3,  2,  if 


Possible  values 


NV  --  any  positive  integer  ^  Gh    (this 

restriction  can  be  easily  removed'* 

NLEV  --  a  list  of  NV  positive  integers 


These  parameters  specify  the  number  of  variables  (NV) 
and  the  number  of  levels  NLEV(i)  for  each  variable  xi  which 
is  used  in  the  description  of  the  input  events.   Since  the 
smallest  variable  value  allowed  is  0,  NLEV(i)  (the  number 
of  levels  for  the  i™  variable)  should  be  one  larger  than 
the  largest  possible  value  of  the  i^n  variable.   In  the 
example  above,  there  are  k   variables:  variable  xl  has  5 
possible  values  (i.e.,  the  domain  of  xl  consists  of  0,1,2,3, 
and  h) ;  x2,  x3,  and  xh   have  3,  2  and  k   possible  values,  respectively. 


•  NCL  NE(NCL) 
Example : 


3  5,  1,  2 


Possible  values:   NCL  --  any  positive  integer 

NE  --  a  list  of  NCL  positive  integers 
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These  parameters  specify  the  number  of  event  sets  (NCL) 
which  are  in  the  input  stream  and  the  number  of  events  for 
each  set  (NE) .   In  the  above  example,  there  are  3  event  sets, 
the  first  has  5  events,  the  second  has  only  1  event,  and  the 
third  has  2  events . 


•  OLIST  (NCL) 
Example : 


10  2 


Possible  values:   some  ordering  of  the  integers  0, 1, . . .,  NCL-1 

(NCL  as  specified  above) 

This  parameter  specifies  the  way  in  which  the  input 
event  sets  are  assigned  symbols  pi,  i=0, 1,  .  ..,  NCL-1 .   The 
consecutive  events  sets  in  the  input  stream  are  numbered 
0, 1,2, ..., NCL-1.   The  first  number  in  the  OLIST  tells  which 
event  set  in  the  input  stream  is  assigned  symbol  f^CL-1^  the 
ith  number  tells  which  event  set  is  assigned  symbol  FNCL-i . 
In  the  example,  the  second  (numbered  as  l^i  event  set  is 
assigned  symbol  F2,  the  first  (numbered  as  0)    is  assigned 
symbol  F  ,  and  the  third  (numbered  as  2")  is  assigned  F^. 


NCRIT  CLIST  TLIST 
Example : 


3  1,  2,  if  0,  0,  0.1 


Possible  values:   NCRIT  --  a  positive  integer  between  1  and  7 

CLIST  --  a  list  of  NCRIT  positive  integers 
between  1  and  7 

TLIST  --  a  list  of  NCRIT  real  values  from  the 
interval  [0,1] 

These  parameters  specify  the  optimality  functional^ 12 
which  is  to  be  used  in  finding  the  'best'  complex  Lq.   NCRIT 
is  the  number  of  cost  criteria  which  are  used.   CLIST  specifies 
criteria  to  be  used  by  listing  the  criteria  numbers  in  order  of 
criteria  preference.   The  criteria  numbers  are  numbers  assigned 
to  each  criterion  as  described  below.   TLIST  specifies  the 
tolerance  which  is  to  be  allowed  for  each  respective  cost  criterion. 

The  cost  functional  is  used  to  select  the  'best'  complex 
from  a  star  (a  set  of  complexes  which  are  maximal  under  inclusion 
and  which  cover  the  same  event  of  Fl  but  no  events  of  F°) .   The 
process  proceeds  as  follows:  each  criterion  is  selected  in  turn 
beginning  with  CLIST(l) .  Using  the  selected  criterion,  a  value 
is  assigned  to  each  complex. 
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in  the  star.  The  minimum  (MINVAL)  and  the  maximum  (MAZVAL) 
cost  of  each  complex  in  the  star  is  computed.  Using  the 
associated  tolerance  limit  TLIST(i),  a  threshold  UBOUND  is 
calculated. 

UBOUND  =  MINVAL  +  TLIST(  i)*(MAXVAL-MINVAL) 

All  complexes  whose  cost  is  greater  than  UBOUND  are  eliminated 
from  further  consideration.   If  this  leaves  only  one  complex, 
the  process  terminates  and  this  is  the  required  'best'  complex. 
Otherwise,  the  next  specified  criterion  is  selected  and  the 
process  is  repeated.   If  all  specified  criteria  have  been  used 
and  more  than  one  complex  remains,  a  complex  with  minimum 
cost  with  respect  to  the  last  criterion  is  selected  as  the 
'best'  complex  (i.e.,  quasi -extremal ) . 

The  user  specified  criteria  characterize  the  formula 
which  he  wants  to  obtain.   Since  the  program  builds  a  formula 
from  selected  complexes,  instead  of  'user  criteria'  (referring 
to  the  whole  formula)  the  program  uses  'program  criteria'  which 
refer  to  the  individual  complexes.   These  program  criteria 
are  designed  to  best  approximate  the  original  criteria 
(which  refer  to  the  whole  formula) .   For  example,  if 
the  user  specifies  'minimum  number  of  terms'  as  a  criterion, 
then  this  criterion  is  substituted  in  the  program  by  the  criterion 
'maximum  number  of  events  which  a  given  complex  covers  in  the 
remaining  set  of  events  to  be  covered'. 


Below  are  described  currently  available  user  criteria^ 
and  corresponding  to  them  program  criteria  (in  parentheses) : 

Criterion 

#      Notation  Description  of  Attributes  (Criteria) 

1         t  --   the  number  of  terms  in  a  formula 

(The  negative  of  the  number  of  events  in  the  event  set 
which  are  covered  by  the  given  complex  but  not  by  any 
of  the  previously  determined  complexes  (stored  in 
NMQ) .   The  negative  is  used  to  have  the  minimum  value 
for  the  criterion  when  the  complex  covers  the  maximum 
number  of  events.) 


s  --  the  number  of  selectors  in  a  formula 

(The  number  of  selectors  in  a  complex.) 


z  —       sum  of  the  costs  of  the  variables  in  a  formula  (z(v)) 
(The  sum  of  the  costs  of  the  variables  in  a  complex: 

z«  =  £  z(i) 
iel 
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z(i)    __  cost  of  the  variable  x^  which  is   specified 

by  the  Z  parameter  (see  Z  parameter  specification 
below) 

I  —  set  of  indices  of  variables  which  appear  in  the 
complex . ) 

g  --  The  degree  of  generalization 

(The  degree  of  generalization  of  a  complex: 

c(E) 

c(e  n  F) 

c(e)  --  number  of  events  in  the  given  complex  E 

c(E  n  F)  --  number  of  events  covered  by  the  complex  E 

which  are  members  of  the  event  set  F1  which 
is  being  covered.) 

W  --  'Weight  distribution' 

(The  negative  of  the  sum  of  weights  of  events  of  the  event 
set  which  are  covered  by  this  complex  but  not  by  any 
previous  complex.   If  this  criterion  is  selected, 
weights  W(K,J)  are  specified  (see  W  parameter  specifi- 
cation below) .   If  the  complex  is  created  to  cover 
event  set  K,    then  the  cost  is  given  by 

-  Z  W(K,J) 
JeE 

E  --  set  of  events  of  the  event  set  k  which  are 

covered  by  this  complex  and  not  by  any  previous 
complex. ) 

I     __  Total  length  of  references 

(The  length  of  references  in  a  complex.   This  is 
the  sum  of  the  number  of  constants  which  appear  in  the 
reference  of  each  of  the  selectors  of  the  complex 
(in  expanded  form) 

Z  i(xi) 
iel 

I  --  as  in  3  above 

l(x±)    --  number  of  constants  in  the  reference  of 
the  selector  for  the  variable  xi .  ) 


The  relative  scope  of  references 

(The  relative  scope  of  references  in  a  complex.   This 
is  the  sum  of  the  mean  deviations  for  variables  in 
the  complex.   The  mean  deviation  of  a  variable  is  the 
average  absolute  difference  between  the  reference 
constants  and  the  mean  of  these  constants.   That  is, 

cost  =  Z   (MD. ) 
lel 


11+ 


MDt  =  Z 


|avgi  -  r 


TCT 


I  --  set  of  indices  of  variables  which  appear  in  the 
complex 

avgj_  --  average  value  of  the  constants  in  the  reference 
for  variable  xi  in  the  complex 


Rj_  --  the  set  of  constants  r  in  the  reference  of  the 
variable  x^  in  the  complex 

c(Rj_)  --  the  number  of  elements  of  R^ 


0   0   1 

0  2   0 

13   2 

103 

113 

2   0  2 

Example : 

Possible  values:   a  set  of  nonnegative  integers  which  properly 

represent  each  event 

F  is  a  list  of  event  sets  F1  defining  the  given  VL  function(s)  . 
F  must  be  completely  specified  (no  default  values).   If  INFORM  = 
'VECTOR',  then  each  event  must  be  represented  as  a  set  of  NV  (number 
of  variables)  integers  each  of  which  specifies  the  respective 
variable  value  for  that  event.   The  value  for  variable  xi  must 
be  less  than  NLEV(i) .   If  INFORM  =  'GAMMA',  then  the  gamma 
representation  for  each  event  must  be  given. 

The  example  above  could  be  the  representation  for  6  events 
in  'VECTOR'  format  with  NV  =  3- 

The  following  parameters  (Z  or  W)  are  specified 
only  if  the  cost  criterion  3  or  5  are  specified, 
respectively.   They  are  specified  in  PL/l  data  format 
and  each  specification  is  terminated  with  a  semi-colon  (;). 


Example:    Z(l)  =  2   Z(k)   =1.2; 

Possible  values:   any  real  numbers 

This  parameter  is  specified  only  if  criterion  3  is 
specified  and  specifies  the  cost  of  each  variable  (i.e., 
the  cost  of  xi  is  specified  by  Z(l)).  Any  variables  for 
which  no  specification  is  given  are  assigned  cost  =  1. 
In  the  example  given,  variable  xl  has  cost  2,  variable  xk 
has  cost  1.2,  and  the  other  variables  have  cost  =  1. 
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•  w 

Example:    W(0,2)  =1.3   W(l,l)  =  2; 
Possible  values:   any  real  numbers 

This  parameter  is  specified  only  if  criterion  5  is  specified. 
W(K,J)  specifies  the  weight  of  Jth  (j=l,2,5,  ...)  event  in  the 
(K+l)s"t  (R=0,1,  ...)  event  set  in  the  list  F.   Any 
event  for  which  no  value  is  given  is  assigned  weight  =1.   in 
the  example  above,  the  second  event  of  the  first  event  set  is 
assigned  weight  =  1.3,  the  first  event  of  the  second  event  set 
is  given  weight  =  2.0,  the  rest  of  the  events  are  given  weight  =  1.0, 
(Note  event  set  numbering  begins  at  0,  event  numbering  within  each 
event  set  begins  with  1.) 

This  concludes  the  parameters  specification  for  one  problem. 
Several  problems  may  be  run  together  by  specifying  all  of  the  parameters 
for  each  problem  in  succession.   No  parameter  values  (except  the  default 
control  parameter  values)  are  carried  over  from  one  problem  to  the  next. 


3-   OUTPUT  PARAMETERS 

Some  of  the  output  parameters  of  the  program  are  reflections  of 
the  input  parameters .   The  following  parameters  are  generated  by  the  program 
(the  text  enclosed  in  the  dotted  line  is  the  actual  information  printed  on 
the  printout  from  the  program) : 

•  Unused  core 

i 1 

Example:   J  AMOUNT  OF  UNUSED  CORE  =  ^k 

This  parameter  specifies  the  amount  of  core  memory 
which  was  estimated  but  not  used. 


•  COMPLEX: 

! : . -  ; . -  ; :  —  — .  -i 

I 


Example:    COMPLEX:   (xl  =  0  2)  (x3  =  l)  (9,  6,  2,  \2\ 


This  parameter  specifies  the  actual  complexes  which 
were  synthesized.  Complexes  of  one  formula  (or  in  the  case 
of  'VL'  mode,  complexes  which  have  the  same  coefficient) 
are  grouped  together  under  the  heading  identifying  the  event 
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class  which  they  cover.   Together  with  each  complex  is 
printed  a  list  of  h   numbers  (COV,  NEW,  IND,  TOT) . 


where : 


COV  --  denotes  the  number  of  events  covered  by  the 
complex  in  the  given  event  set  F1  (i.e., 
event  set  associated  with  the  current  decision 
class) 

NEW  --  is  equal  to  COV  minus  number  of  events  in 
F^  covered  by  the  previous  complexes  on  the 
output  list  ('new  covered'). 

IND  --  is  equal  to  COV  minus  number  of  events  in 

F^  covered  by  all  the  other  complexes  in  the 
cover  ('independently  covered') 

TOT  —  is  equal  to  the  total  number  of  events  in  F^ 

which  are  covered  by  this  complex  and  previous 
complexes  in  the  cover 


•  DELTA 

,—  —  —  —  —  —  —  —  —  —  _  —  —  —  —  — 

Example:   [  DELTA  for  this  set  is  0 


This  parameter  specifies  an  estimate  (upper  bound)  on 
the  distance  (measured  in  number  of  complexes)  of  the  formula 
produced  from  the  minimal  formula.   In  the  above  example, 
since  DELTA  is  0  the  set  of  complexes  is  minimal.   Note  that 
if  stars  are  trimmed  (i.e.,  an  intermediate  star  exceeds 
MAXSTAR  complexes)  then  DELTA  is  no  longer  a  valid  estimate. 
If  complexes  are  a  result  of  trimmed  star,  a  message  to  this 
effect  is  printed. 


INTERMEDIATE  STAR  SIZE 

r—— —  — ___ —  _ 

Example:   !  Largest  intermediate  star  size  is  5& 


This  parameter  specifies  an  estimate  of  the  maximum 
number  of  complexes  which  are  stored  in  the  star  list.   The 
NGE  parameter  may  be  reduced  to  approximately  MAXSTAR  plus 
this  estimate  if  this  problem  is  run  again. 


k.      EXAMPLES  OF  PROGRAM  INPUT  AND  OUTPUT 

Two  examples  of  the  input  and  output  formats  for  the  program 
AQVAL/i  (AQ7)  are  presented  below.   The  first  example  is  further  illustrated 
with  GLD  representations  of  the  formulas  produced  by  the  program.   Included 
with  each  line  of  the  input  examples  is  a  description  which  lists  the  names 
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of  the  parameters  which  are  specified  on  that  line,  a  brief  indication 
of  the  significance  of  each  parameter  and  the  actual  values  which  the 
parameter  assumes  in  the  example. 

The  first  few  lines  of  each  example  of  program  input  is  the 
JCL  specification.   These  cards  must  be  included  before  any  parameter 
specification  to  gain  access  to  the  program  and  distribute  the  program 
output  properly.   Preceding  the  JCL  specification  must  be  the  ID  cards 
which  contain  the  time  and  region  estimates.  The  amount  of  these 
resources  which  the  program  requires  is  a  complex  function  of  the 
control  parameters  and  the  dimension  and  size  of  the  event  sets  which 
are  to  be  covered.   Samples  of  various  parameters  and  the  amount  of  time  and 
region  to  successfully  complete  the  execution  of  the  program  are  given  in 
Figure  1 .   In  general,  for  small  problems,  the  region  estimate  must  be 
at  least  ll6  k  and  the  time  estimate  can  be  less  than  30  seconds.   For 
very  large  problems,  a  region  estimate  of  l80k  and  time  of  6  minutes 
is  probably  sufficient.   Note  that  the  region  estimate  must  appear  on  both 
the  ID  card  and  the  EXEC  card  of  the  JCL.  Additional  blank  lines  may  be 
inserted  for  clarity  in  several  points  in  the  input.  They  must  not  appear, 
however,  between  any  JCL  cards  or  between  the  semi-colon  which  terminates 
the  control  parameters  and  the  title  (unless  the  control  parameter  TITLE 
has  allowed  for  this) . 

The  output  of  the  program  begins  with  the  title  which  was 
specified  in  the  input.   The  values  of  all  control  parameters,  the  type 
of  each  variable  and  event  set  sizes  are  listed.   Then,  the  amount  of 
core  remaining  is  given,  which  is  the  amount  of  core  which  was  estimated 
but  not  used.   Next,  the  OLIST  and  optimization  criteria  are  listed 
followed  by  the  actual  event  sets  which  were  specified.   The  VL  formula(s) 
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follow.  Each  complex  (or  term)  in  the  synthesized  formula (s)  is 
printed  along  the  DELTA  for  this  set  of  complexes.  DELTA  gives  an 
estimate  (upper  bound)  of  the  number  of  complexes  by  which  this  formula 
exceeds  the  minimal  formula.   If  the  formulas  are  a  result  of  trimmed 
stars  (the  size  of  some  intermediate  star  exceeded  the  MAXSTAR  parameter), 
a  message  is  printed  which  states  that  the  complexes  are  a  result  of 
trimmed  stars  and  thus  the  value  of  DELTA  is  only  an  estimate  (not  an 
upper  bound)  of  the  distance  of  the  formula  generated  to  the  minimal 
formula.   The  largest  star  size  and  the  largest  intermediate  star  size 
are  printed.   The  largest  intermediate  star  size  plus  the  value  of  MAXSTAR 
gives  an  estimate  of  the  minimum  size  of  the  parameter  NGE  which  was 
actually  necessary  in  the  synthesis  of  this  set  of  complexes.  After  the 
entire  set  of  formula(s)  has  been  printed,  a  message  indicating  a  normal 
termination  is  given  which  indicates  the  successful  completion  of  the 
synthesis  of  formulas  for  this  problem. 

The  first  example  demonstrates  the  three  modes  of  operation  of 
the  program  as  specified  by  the  parameter  MODE.   In  addition  to  the  input 
stream  and  the  output  generated,  the  formulas  which  were  generated  are 
illustrated  in  Fig.  2,  3  and  k   by  Generalized  Diagrams  (GLD) .   The  second 
example  demonstrates  the  use  of  criteria  3  and  5,  the  specification  of 
the  associated  cost  and  weight  parameters  (Z  and  w),  and  the  different 
tolerance  specifications  for  each  criterion. 

GLD  Description  for  Example  #1 

The  GLD  is  a  graphical  model  of  the  event  space  and  the 
complexes  which  cover  the  events  in  this  space.   Each  square  on  the  GLD 
represents  one  event  in  the  event  space;  a  set  of  events  enclosed  in  an 
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area  or  set  of  areas  can  represent  a  complex.   In  this  GLD,  the 
variables  xl  and  x2  specify  the  row  in  which  a  particular  event  is  to 
be  placed,  the  variables  x3  and  xk   specify  the  column.  An  event  set 
is  a  set  of  squares  with  the  same  number.   If  the  event  sets  are  to  be 
covered  in  VL  mode,  then  the  numbers  (i)  in  the  squares  correspond  to 
the  set  F  and  to  the  terms  of  the  DVL  formula  with  coefficient  i.   If  the 
event  sets  are  to  be  covered  in  IC  or  DC  mode,  then  the  numbers  (i) 
correspond  to  those  events  with  a  function  value  of  1  for  the  i^n 
function  and  a  function  value  of  0  for  the  other  functions.   The  events 
which  are  given  to  the  program  are  indicated  in  the  diagrams  by  numbers 
in  the  squares.   The  complexes  of  each  type  of  cover  (VL,  IC,  DC)  which  were 
generated  by  the  program  are  illustrated  in  the  respective  diagrams 
(Fig.  2,  3,  k)   by  open  or  shaded  areas  together  with  their  VL-]_ 
representation.  Some  observations  about  these  three  different  covers 
follow . 

VL  mode:   The  complex  which  covers  the  event  set  F^  does  not 

cover  events  of  any  of  the  sets  F1,  1=3,2,1,0.   The  complex  which  covers  the 

2  h.  x 

event  set  F  ,  however,  also  covers  some  events  of  sets  F  and  FJ .   Since  this 

formula  is  evaluated  by  testing  the  complexes  associated  with  event  sets 

F^"  and  F^  before  the  complex  associated  with  the  set  F  ,  the  events  in 

the  original  event  sets  are  classified  correctly. 

IC  mode:   Each  complex  does  not  cover  any  events  from  any  event  sets 

other  than  the  one  which  it  was  specified  to  cover.   The  complexes  which 

cover  event  sets  F  and  W   intersect  with  the  complex  which  covers  event 

set  F  in  an  unspecified  area  of  the  event  space.   It  is  thus  possible 

that  an  event  in  this  area  could  satisfy  two  different  formulas.   (Such 

a  situation  could  be  interpreted  as  an  uncertain  classification.) 
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DC  mode:  The  complexes  are  similar  to  those  produced  in  IC 
mode.  However,  they  do  not  intersect  in  any  area  of  the  space  hence  at 
most  one  formula  is  satisfied  by  the  set  of  complexes  for  any  event  in 
the  event  space. 
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5-   PROGRAM  IMPLEMENTATION 
5-1  General  remarks 

The  basic  algorithm  used  in  the  program  for  the  synthesis  of 
DVL-l  formulas  is  algorithm  Aq  described  in  papers  '    '   .    In  this 
section,  we  will  describe  all  the  basic  concepts,  data  structures  and 
the  function  of  each  procedure  of  the  program.  This  description  will  be 
sufficient  for  general  understanding  of  the  program  and  also  can  serve  as 
a  guide  for  the  detailed  study  of  the  program  implementation.   The  complete 
listing  of  the  program  is  given  in  Appendix  B  of  this  paper. 

5.2  Program  representation  of  complexes  and  events 

All  events  and  complexes  are  represented  in  the  program  by 

PL/l  bit  strings.  Each  variable  x.  in  an  event  or  complex  is  assigned 

a  different  substring  of  this  bit  string.  Each  position  of  the  substring 

corresponds  to  a  value  of  this  variable.  Consequently,  the  length  of  the 

substring  assigned  to  variable  x.  equals  the  number  of  values  (NLEV(i)) 

in  the  domain  of  this  variable  and  the  length  of  the  entire  string  is 

the  sum  of  the  lengths  of  each  substring,  i.e. 

NV 

Z  NLEV(i) 
i=l 

(NV  --  number  of  variables) 

The  leftmost  position  of  the  substring  corresponds  to  value  0  and  the 
rightmost  to  value  d\,  i.e.,  maximum  value  of  the  given  variable  xi .   If 
in  an  event  (complex)  a  variable  x±   has  (is  referred  to)  value  j,  then 
(j+l)th  position  of  the  substring  is  set  to  value  1.   Remaining  positions 
are  set  to  value  0.   For  example,  the  event 

e  =  (0,  3,  2,  h) 


ko 


from  event  space  E(5,  lh    ^,    5)  (i.e.,  NLEV(l)  =  3,  NLEV(2)  =  k, 
NLEV(3)  =  k,    NLEV(4)  =  5)  is  represented  as: 


'100  0001  0010  00001' B 
xl    x2   x3    x*4- 


CTearly,  each  substring  in  a  string  representing  an  event 
has  exactly  one  bit  set  to  1,  and  each  substring  in  a  string  representing 
a  non-empty  complex  has  at  least  one  bit  set  to  1.   If  in  a  complex  a 
given  variable  is  absent,  then  the  substring  corresponding  to  this  variable 
has  all  bits  set  to  1.  For  example,  the  complex 

(xl=2)(x4=l,3) 


is  represented  as; 


'001  1111  1111  01010': 
xl   x2    xj     xh 


5 .3   Internal  structure  for  storing  intermediate  stars 

A  star  G(e|E)  of  an  event  e  against  event  set  E  is  generated  in 
steps,  each  step  corresponds  to  an  event  in  E  (more  precisely,  each  step 
involves  the  multiplication  of  the  current  set  of  complexes,  called  an 
'intermediate  star',  by  complexes  resulting  from  the  extension  of  e  against 
the  consecutive  event  e^  from  E) .   Thus,  an  intermediate  star  IG(e |e)  is  a 
set  of  maximal  complexes  which  cover  e  and  do  not  cover  some  of  the  events 
from  E  (rather  than  all  events,  as  would  be  the  case  with  G(e|E)). 

An  intermediate  star  is  represented  by  a  doubly  linked  list, 
each  element  of  which  contains  a  pointer  to  the  previous  and  to  the 
next  element  in  the  list.   In  addition  to  a  forward  pointer  (FPTR)  and 


a  backward  pointer  (BPTR)  each  element  of  the  list  contains  the  bit 
string  representation  of  the  complex  (num),  and  a  field  (VAL)  reserved 
for  the  value  of  any  criterion  used  to  evaluate  this  complex  during  the 
selection  of  quasi-extremal  (i.e.,  selection  of  'best  complex'). 

The  list  can  be  accessed  by  a  pointer  which  points  to  the  first 
element  of  the  list  (the  head  of  the  list)  and  a  pointer  which  points  to 
the  last  element  of  the  list  (the  tail  of  the  list) .   In  the  process  of 
star  generation  three  (doubly  linked)  lists  are  used: 

1.  a  list  representing  the  last  intermediate  star  (with 
head  FGPTR  and  tail  BGPTR), 

2 .  a  list  representing  the  intermediate  star  under 
construction  (with  head  and  tail  pointers  NFGPTR 
and  NLGPTR), 

3.  a  list  of  free  elements  (with  head  and  tail  pointers 
NEXTFG  and  LASTFG),  reserved  for  additional  complexes. 

All  3  lists  are  stored  in  an  array  structure  G.  A  graphical 

representation  of  storing  an  intermediate  star  within  the  structure  G 

follows: 


FGPTR 


118 


G(ll8) 


G(113) 


G(115) 


BGPTR 


110 


G(110)   ,. 
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complex  1 

complex  2 

complex  3 

complex  k 

criterion 
value  1 

criterion 
value  2 

criterion 
value  3 

criterion 
value  k 
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5 .k     Structure  of  the  program 

The  first  set  of  PL/1  statements  (approximately  150  statements) 
allocate  the  proper  storage  for  all  global  variables.   The  input  parameters 
are  read  and  their  values  echoed  on  the  printout.   The  type  of  each  variable 
is  recorded  in  a  bitstring  (TYPCOV)  where  'O'B  in  the  ith  position  indicates  that 
variable  Xj_  is  of  type  INTERVAL,  and  '1'  that  it  is  of  type  FACTOR.   The 
events  are  converted  to  bitstring  representations  as  described  above  and 
loaded  into  an  array  E  so  that  each  row  of  E  corresponds  to  an  event  set, 
and  each  column  corresponds  to  an  event  in  each  set. 

The  next  set  of  statements  (approximately  40  statements)  controls  the 
selection  of  event  sets  for  covering.   The  selection  of  these  sets  depends  on 
the  mode  of  operation  of  the  program  (specified  by  the  parameter  MODE)  and 
on  the  order  in  which  the  sets  are  to  be  covered  (as  specified  by  the 
parameter  OLIST) .  A  row  of  the  array  E  corresponding  to  the  set  of  events 
to  be  covered  is  selected  and  the  array  FDOWN  is  loaded  with  the  set(s)  of 
events  (and  complexes,  in  the  case  of  DC  mode)  against  which  the  covers 
are  to  be  generated.   In  addition,  two  bit  arrays  F1N0MQ  and  FBfOSTAR  are 
maintained.   These  arrays  are  in  a  one-to-one  correspondence  with  the 
events  which  are  to  be  covered  (the  columns  in  the  selected  row  of  E), 
(the  array  F1N0MQ,  indicates  elements  which  belong  to  the  set  E-,  and  the 
array  F1N0STAR  indicates  the  elements  which  belong  to  the  set  E^  in  the 
algorithm  A  presented  in  Figure  1  of  paper  ) .  A  one  ('l'B)  in  a  position 
of  F1N0MQ,  indicates  that  the  corresponding  event  of  E  belongs  to  E-, , 
i.e.,  it  has  not  yet  been  covered  by  any  quasi-optimal  complex  (L^)  ;  a 
one  ('l'B)  in  a  position  of  F1N0STAR  indicates  that  the  corresponding 
event  belongs  to  E  ,  i.e.,  it  has  not  yet  been  covered  by  any  star  (i.e., 
by  any  of  the  complexes  of  previously  generated  stars) . 
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Next  in  the  program  are  the  loops  controlling  the  selection 
of  the  events  to  he  covered  according  to  the  arrays  F1N0MQ,  and  F1N0STAR 
(approximately  50  statements) .   In  the  first  loop,  an  event  which  has 
not  been  covered  by  any  previous  star  (i.e.,  events  of  Ep)  is  selected 
according  to  FINOSTAR.  A  star  is  generated  for  this  event  and  the 
'best  complex'  (quasi-extremal  Lq)  is  extracted  by  calls  to  the  procedures 
STAR  and  BESTLQ,.   The  resulting  qua  si -extremals  are  stored  in  a  list  MQ.INT 
along  with  a  bit  MQ.TRIM  which  indicates  whether  this  complex  was  a  result 
of  a  trimmed  star  or  not.  Once  the  set  of  events  represented  by  FINOSTAR 
(i.e.,  events  of  Ep)  has  been  exhausted,  the  process  is  repeated  with  the 
set  represented  by  F1N0MQ  (i.e.,  set  Ex) .   The  number  of  passes  through 
this  second  loop  is  calculated  and  stored  in  the  parameter  DELTA  as  an 
estimate  of  the  distance  of  this  set  of  formulas  from  the  optimal  set  (see 
equation  (12)  in  paper2).   The  resulting  set  of  complexes  is  printed  by 
a  call  to  the  procedure  PRCOVER.   Next  a  new  set  of  events  to  be  covered 
(row  of  E)  is  selected  and  the  array  FDOWN  is  loaded  with  new  sets  of 
events  (and  complexes  in  DC  mode)  against  which  a  new  cover  is  to  be 
generated.  After  all  event  sets  which  were  to  be  covered  have  been 
covered  and  the  complexes  printed  out,  the  next  set  of  parameters  are 
read  in.   If  no  more  parameters  exist,  the  execution  is  stopped. 

5-5  Basic  procedures  in  the  program 

The  remainder  of  the  program  consists  of  a  set  of  procedures 

each  of  which  performs  a  specific  task  as  is  explained  below: 

STAR  —  This  is  the  procedure  which  computes  a  star 

G(e|FD0WN)  of  an  event  e  selected  from  the  given 

row  of  the  array  E  (representing  event  set  to  be 

covered)  against  the  event  set  (and  set  of 

complexes  in  the  DC  mode)  represented  by  FDOWN. 

The  star  is  generated  in  steps,  each  step 

corresponding  to  one  element  of  FDOWN  and  producing 

an  intermediate  star.   In  each  step,  first  the  procedure 


kk 


MIX  is  used  which  generates  an  'elementary  star' 
G(e|L),  where  L  is  the  selected  element  of  FDOWN, 
i.e.,  set  of  maximal  complexes  which  cover  e  and  do  not 
cover  L.   Now,  the  last  intermediate  star  is  multiplied 
by  the  'elementary  star'  G(e|L)  and  absorption  rules 
are  used  before  yielding  the  new  intermediate  star. 
The  process  ends  when  the  last  element  of  FDOWN  has  been 
used.   If  any  of  the  intermediate  stars  have  more 
elements  than  MAXSTAR,  a  procedure  TRIMSTAR  is  called. 


MIX  —  This  procedure  produces  the  star  G(e|L)  where  L  is  an 

element  (event  or  complex)  in  FDOWN.   The  stars  produced 
by  MIX  are  represented  by  the  arrays  INTP  and  POSP.  Each 
complex  of  a  star  of  an  event  against  another  event  or 
complex  involves  only  one  variable  (i.e.,  has  only  one 
literal),  therefore  in  this  type  of  star,  only  the  values 
of  one  variable  must  be  given  in  the  representation  of  a 
complex  in  the  star.   In  the  program,  INTP(i)  is  a 
bitstring  which  represents  the  values  of  the  variable 
whose  index  is  stored  in  POSP(i)  (only  non-empty  complexes 
are  retained)  and  NUM  indicates  the  number  of  non-empty 
complexes  which  were  placed  into  the  star.   The  format  of 
each  of  these  complexes  produced  by  MIX  depends  on  the 
type  of  the  variable  which  generates  the  complex.   The 
procedure  is  as  follows:   the  complement  of  the  literal 
of  the  event  or  complex  of  FDOWN  against  which  a  cover 
is  being  generated  is  computed.   If  this  does  not  contain 
the  value  which  the  variable  assumes  in  the  chosen  event  e 
of  E  then  the  complex  is  thrown  out  (i.e.,  because  it  is 
empty) .   If  the  type  of  this  variable  is  FACTOR,  the 
complemented  literal  is  used  as  the  literal  of  the  complex 
and  loaded  into  the  next  place  of  INTP.   Otherwise,  the 
interval  of  this  complement  which  contains  the  value  of  the 
corresponding  literal  of  the  chosen  event  e  is  placed  into 
INTP.   In  each  case  where  there  is  a  non-empty  complex,  the 
index  of  the  variable,  which  corresponds  to  the  literal 
in  INTP,  is  placed  into  POSP. 


TRIMSTAR  --  This  procedure  is  invoked  by  the  procedure  STAR  when 
the  size  of  an  intermediate  star  becomes  greater  than 
MAXSTAR.   It  reduces  the  star  to  the  CUTSTAR  best  complexes 
according  to  criteria  1  and  2  defined  above  under  the 
parameter  CLIST.   (These  criteria  measure  the  number  of 
events  in  the  set  to  be  covered  (indicated  by  F1N0MQ)  and 
number  of  literals  in  a  complex,  respectively) . 


NUMLIT  --  This  procedure  counts  the  number  of  literals  in  a  given 
complex . 
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LCOVER  --  This  procedure  counts  the  number  of  events  covered 

by  a  given  complex  and  not  covered  by  any  previously  chosen 
complex  (i/l)  (i.e.,  the  number  of  events  indicated  by 
F1N0MQ,  which  are  covered  by  the  complex)  . 


BESTLQ  --  This  procedure  uses  the  criteria  specified  in  the  input 
to  select  the  'best'  complex  (quasi-extremal  Lq)  from 
the  given  star.   The  selection  process  is  described  in 
section  2  (parameter  CLIST) . 


CRITVAL  --  This  is  a  function  which  returns  the  value  of  a 

criterion  applied  to  a  given  complex.   The  value  which 
is  returned  by  CRITVAL  for  each  criterion  is  described  in 
section  2  (parameter  CLIST) . 


RANGE  --  is  a  procedure  internal  to  CRITVAL  used  to  determine  the 
number  of  values  of  a  variable  in  the  given  complex.   0.0 
is  returned  if  the  variable  does  not  appear  in  the  complex 
(i.e.,  the  variable  accepts  all  possible  values). 

VARIANCE  --  is  a  procedure  internal  to  CRITVAL  used  to  determine 
the  mean  deviation  (described  in  criterion  7  in  section  2 
under  the  parameter  CLIST) of  the  values  of  a  variable  in 
the  given  complex.   0.0  is  returned  if  the  variable  does 
not  appear  in  the  complex. 


PRTSTAR  --  prints  out  the  elements  of  a  star, 
PRTQLQ  --  prints  out  trace  information  on  L<1. 


READVEC  --  reads  input  events  in  vector  format  and  converts  them 
to  the  internal  bit  string  representation. 


READGAM  --  reads  input  events  in  gamma  format  (i.e.,  as  event 
numbers)  and  converts  them  to  internal  bit  string 
representations . 


CHECKE  --  is  an  auxilliary  procedure  used  by  STAR  in  the  process 
of  combining  stars  (application  of  absorption  laws) . 


CKBE  --  is  an  auxilliary  procedure  used  by  STAR  in  the  process 
of  combining  stars. 
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REMOVEGN(REMOVEG)  --  is  used  to  remove  a  complex  from  the 

star  during  the  generation  and  selection  processes. 


GETG  —  obtains  an  element  of  G  from  the  free  list  for  use  in  a 
star. 


RETURNG  --  returns  a  list  of  elements  of  G  to  the  free  list. 


EINCPLX  --  is  a  function  which  determines  if  a  given  event  in  the 
set  to  be  covered  is  covered  by  a  given  complex. 


PRCOVER  --  prints  out  a  cover  in  a  readable  format. 


FRECORE  --is  used  by  the  initialization  procedure  to  determine 
the  amount  of  available  core  remaining.   Note  that  it  is 
dependent  on  the  method  of  implementing  pointer  variables 
and  on  the  layout  of  OS/360  control  blocks. 


5.6  Remarks  about  other  related  programs 

The  program  AQVAL/l  (AQ7)  is  currently  the  main  and  most 
developed  program  related  to  the  computer  implementation  of  VL-.  .  Other 
programs  which  are  currently  available  include: 

•  OMNIBUS  --  A  program  which  can: 

A.  Evaluate  a  given  VL-j_  formula  (in  which  complexes 
are  represented  as  bit  strings)  for  any  given 
event . 

B.  Quantize  a  given  set  of  variables  according  to 
a  set  of  given  thresholds. 

C.  Select  from  the  original  set  of  events,  subsets 
of  events  such  that 

(i)  they  can  involve  only  certain  variables 
(out  of  all  available  variables) 

(ii)  they  can  constitute  any  contiguous  subset 
of  the  original  set  of  events. 


(iii)  gather  statistics  of  performance  of  VL-i 
formulas  used  as  classification  rules*. 


/ 


This  last  function  is  implemented  in  an  extended  form  in  the  program 
CONFUS  (see  next) . 
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UNICLASS  (AQ8) 


A  special  program  developed  for  the  synthesis 
of  VLj  formulas  from  event  sets  which  is  able 
to  create  exact  or  approximate  (with  a  user 
specified  degree  of  approximation)  covers  of 
a  set  of  events  against  its  complement . 21 


•     AQVAL/l    (AQ9) 


A  program  which   synthesizes   quasi-optimal  DVL-. 
formulas   from  non-optimal  DVLi   formulas . 


CONFUS 


A  special  program  implementing  various  methods 
of  evaluation  of  families  of  DVLi  formulas  and 
printing  a  confusion  matrix  which  relates 
evaluated  values  to  the  correct  known  values 
for  the  formulas  (appl:  for  testing  formulas  used 
as  decision  rules  for  pattern  recognition  problems) 


•  HISTOGRAM 


A  program  which  prints  a  graph  which  illustrates 
the  distribution  of  frequencies  of  individual 
values  of  a  variable  in  an  event  set.  (appl; 
determination  of  thresholds  for  quantization 
units  in  designing  learning  and  testing  for 
pattern  recognition  problems) 


•  SIM1 


A  program  which  detects  symmetry  in  (incompletely 

specified)  VL]_  functions . 

(This  program  is  still  under  development.) 
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